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ABSTRACT 

Phase-resolved high resolution optical spectroscopy has revealed narrow NIII and 
Hell emission lines from the soft X-ray transient Aquila X-l during its 2004 outburst 
that move as a function of the orbit consistent with the phasing of the donor star. 
Under the assumption that these lines come from the irradiated side of the donor star, 
we can constrain its K% velocity to ^247±8 km s _1 , and derive a mass function of 
/(A/i)^1.23±0.12 M@. Estimates for the rotational broadening based on the emission 
components suggest a possible massive neutron star of ^1.6M Q (at 95% confidence). 
However, an updated ephemeris and additional high resolution spectroscopy of AqlX-1 
during a future outburst are warranted in order to confirm that the narrow lines indeed 
originate on the donor star surface, and reliably characterise the system parameters of 
this important X-ray binary. Spectra taken during the end of the outburst show that 
the morphology of the emission lines changed dramatically. No donor star signature 
was present anymore, while the presence of narrow low-velocity emission lines became 
clear, making AqlX-1 a member of the slowly growing class of low- velocity emission 
line sources. 

Key words: accretion, accretion disks - stars:individual (Aquila X-l) - X- 
rays:binaries. 



1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are systems where a com- 
pact object, either a neutron star or black hole, accretes 
material via Roche-lobe overflow from a low mass (^IMq) 
companion. Soft X-ray transients (SXTs) are an important 
sub-class of the LMXBs. Contrary to the persistent LMXBs, 
that always accrete material at high rates, they occasionally 
have short outbursts (~weeks up to months) before return- 
ing back to quiescence. These systems are especially interest- 
ing for the prospect of quiescent studies of the mass donor, 
since they provide the only direct method to obtain orbital 
solutions and system masses for LMXBs (e.g. Charles & Coe 
2005). Such system parameter studies are usually not pos- 
sible during an outburst because the photospheric emission 
from the companion star is swamped by the emission from 
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the X-ray irradiated disk. Unfortunately, in quiescence many 
SXT companion stars are too faint (U>22-23) to allow kine- 
matic studies. This has substantially hampered detailed in- 
terpretation of their evolution, and our ability to construct 
compact object mass distributions. 

AqlX-1 is a SXT that shows one of the shortest recur- 
rence times, making it one of the most intensively studied 
systems. AqlX-1 shows quasi-periodic outbursts, with ap- 
proximately one year intervals (Priedhorsky & Terrell 1984). 
During such an outburst it shows Type I X-ray bursts, indi- 
cating that the compact object is a neutron star (Koyama 
et al. 1981; Czerny et al. 1987). Its optical counterpart was 
detected by Thorstensen et al. (1978) during an outburst at 
a magnitude of B~17. An 18.97 hr periodic modulation in 
the optical was observed during an outburst, which was in- 
terpreted as the orbital period (Chevalier & Ilovaisky 1991). 
In quiescence the counterpart fades to U~21.6, but is only 
0.46" away from a U=19.26 star (Chevalier et al. 1999), 
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Figure 1. Top: X-ray light curve of AqlX-1 obtained by the 
All Sky Monitor onboard RXTE showing the 2004 outburst. The 
arrows indicate the times of our optical spectroscopy. Bottom: 
X-ray hardness ratio defined as 3.0-12.1 keV/1. 3-3.0 keV. 

severely hampering kinematic studies of AqlX-1. Despite 
intense studies there are currently no strong constraints on 
the system parameters. 

Photometric observations show that the light curve 
is double humped during quiescence, changing to single 
humped in outburst (e.g. Welsh et al. 2000; Chevalier & Ilo- 
vaisky 1991,1998; Shahbaz et al. 1998). This is interpreted 
as ellipsoidal variations from the Roche-lobe filling donor 
star in quiescence, changing in outburst to reflection due to 
heating of the side of the secondary facing the neutron star. 
This suggests that during an outburst the optical light is 
dominated by the irradiated side of the secondary (Welsh 
et al. 2000). The spectrum of AqlX-1 in quiescence, after 
careful deconvolution of the spectrum of the contamination 
star, shows absorption features that can be associated with 
a K6-M0 V star, as well as emission lines typical for LMXBs 
(Chevalier et al. 1999). 

In this paper we will explore the possibility to determine 
the system parameters of AqlX-1 during an outburst using 
narrow Bowen fluorescence lines. This technique was first ap- 
plied to ScoX-1 (Steeghs & Casares 2002). Narrow emission 
line components coming from the irradiated companion star, 
especially in the Bowen blend (a blend of N III 4634/4640 A 
and CIII 4647/4650 A lines), were used to trace its orbital 
motion and determine the mass function. Since then, this 
technique has led to the discovery of donor star signatures 
in other persistent LMXBs such as X 1822-371, 4U 1636-53 
and 4U 1735-44 (Casares et al. 2003, Casares et al. 2006), or 
transient sources where the optical counterpart is too faint 
in quiescence, such as GX339— 4 (Hynes et al. 2003a). 



2 OBSERVATIONS AND DATA REDUCTION 

From May until June 2004 AqlX-1 was in outburst (see 
Fig.[TJ, and we were able to obtain a total of 28 spectra 
using the FORS-2 spectrograph attached to the VLT/UT4 
(Yepun Telescope) at Paranal Observatory (ESO). In Ta- 
ble[T]we give an overview of all the observations. The spec- 
tra were taken with the 1400V volume-phased holographic 
grism, and a slit width of 0.7", giving a wavelength coverage 



Table 1. Overview of the observations of AqlX-1. Indicated 
are the UT dates of each observation, R magnitude according 
to observations by Maitra & Bailyn (2005), number of observa- 
tions during that night, exposure time, and corresponding orbital 
phases of the observations according to the weighted ephemerides 
by Chevalier et al. (1991,1998) and Welsh et al. (2000). 



Date 


R 


No. 


Exp. 


Orb. phases 


(dd-mm-yy) 


(mag) 




(s) 




26-05-04 


16.3 


1 


900 


0.91 


27-05-04 


16.3 


8 


600 


0.88-0.90,0.14-0.17 


28-05-04 


16.3 


11 


900 


0.17, 0.30-0.43 


18-06-04 


17.5 


8 


674 


0.76-0.84 



of 4514-5815 Awith a resolution of 70 km s _1 (FWHM). The 
seeing during these observations varied between 0.5 and 2.1 
arcsec. During daytime He, Ne, Hg and Cd arc lamp expo- 
sures were taken for the wavelength calibration scale. 

We de-biased and flat-fielded all the images and used 
optimal extraction techniques to maximise the signal-to- 
noise ratio of the extracted spectra (Home 1986). The pixel- 
to-wavelength scale was determined using a 4th-order poly- 
nomial fit to 20 reference lines resulting in a dispersion of 
0.64 Apix -1 and rms scatter <0.05 A. We also corrected for 
velocity drifts due to instrumental flexure (always <5 km 
s _1 ) by cross-correlating the sky spectra using a small re- 
gion around the [O I] sky line at A5577.4. We thus realigned 
all spectra using this single skyline as the zero point. 



3 DATA ANALYSIS 

Fig. [5] shows the average spectrum of both the 26-28 May 
2004 run (bottom) and the 18 June 2004 observations (top). 
We note that the spectra are dominated by high excitation 
emission lines, notably from Hell A4686, Hell A5411, H/3 
and the Bowen region AA4630-4650, as expected from an X- 
ray binary during outburst. Comparing the two spectra in 
Fig-El an d especially the close-ups of the most important 
emission lines in Figs.[3]y4] the first thing we notice is that 
the morphology of the emission lines has changed signifi- 
cantly. During the 26-28 May 2004 run the emission lines 
consist of a broad component (most likely coming from the 
accretion disk), and a narrow component that shows a sig- 
nificant velocity shift between the individual spectra (these 
components are clearest in Fig.[3]|. On the other hand, dur- 
ing the 18 June 2004 observations the lines are exception- 
ally narrow (see Table[2}, and both the broad base and the 
moving component have disappeared. In order to quantify 
the change in emission lines, we have listed the equivalent 
widths, ratio of peak flux over continuum flux and FWHM of 
the main lines derived from average spectra of the May and 
June 2004 observations in Tabled Furthermore, we have 
also indicated the R magnitude, as determined by Maitra 
& Bailyn (2005), in Table[T] The only line that has broad- 
ened a little compared to previous observations is H/3. How- 
ever, we do note that the Hell A4859 emission line from 
the Pickering series is also present (see Fig.|4|, and could 
contribute to this broadening, as was also the case in the 
LMXB XTEJ2123-058 (Hynes et al. 2001). Most impor- 
tantly, the narrow lines in the June data do not show any 
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Figure 2. Average optical spectrum of AqlX-1 during 26 to 28 May 2004 run (bottom), and during the 18 June observations (top). We 
have indicated the most prominent lines. 



orbital motion, despite phase coverage of ~0.1P or b near the 
orbital conjunction. A close inspection of the individual May 
2004 spectra shows that these weak and narrow stationary 
lines are also present but swamped by the much stronger 
disk and moving components. Although these narrow lines 
are obscured most of the time by the moving lines, the data 
is consistent with them being present in all spectra. Unfor- 
tunately, since we do not have an absolute flux calibration 
for the spectra, it is not clear if the line flux from these 
stationary components has evolved between May and June 
2004. 

Taking a closer look at the region around the Bowen 
AA4630-4650 emission in the individual spectra (see Fig.[3j) , 
we notice that it is dominated by narrow emission lines from 
NTH A4640/4634. Curiously, there is no indication of the 
presence of narrow C III lines , as was the case in e.g. Sco X- 
1 (Steeghs & Casares 2002), but we do note that also in 
other LMXBs the N III Bowen lines are stronger than the 
CIII lines (e.g. Casares et al. 2006). Fig.[5]shows that these 
NTH lines are moving, and we interpret this as orbital mo- 
tion of a confined region in the system. Furthermore, a close 
inspection of Hell A4686 shows that there is a strong nar- 
row component (making up ~30% of the total line flux) 
that appears to be moving in phase with the narrow Bowen 
lines (see Fig.[5jl . Given that the off-set between the narrow 
He II and Bowen lines does not change between the spectra 
strengthens our identification of the narrow Bowen lines as 
NIII. 

In order to estimate the orbital phases of our spectra we 
started with the orbital period derived by Chevalier et al. 
(1998) and phase zero by Welsh et al. (2000). Phase zero 
is defined as the minimum in the photometric lightcurve, 
which is expected to correspond to inferior conjunction of 
the secondary if the photometric lightcurve is indeed domi- 
nated by its emission. Unfortunately, this phasing does not 
match up with the published phase zero by Chevalier et al. 
(1991). The differences between these zero points is 0.1 or- 
bital phase, while the error suggests only a difference of 0.04 
orbital phase. Since we cannot be certain what causes this 



Table 2. Overview of the characteristics of the average main 
emission lines observed in AqlX-1 during the May and June 
2004 observations respectively. For each line the equivalent widths 
(EW), ratio of peak (Fp^^) over continuum flux (F con t), and 
FWHM are shown. 



Line 




May 2004 


June 2004 


H/3 


EW (A) 

-^peakZ-^con 

FWHM (km s" 1 ) 


0.52±0.01 
1.2±0.1 
447±6 


1.47±0.04 
1.4±0.1 
561±13 


Hell 


EW (A) 

-^peakZ-^con 

FWHM (km s- 1 ) 


3.49±0.02 
1.3±0.1 

539±5 


2.26±0.06 
1.6±Q.l 
153±6 


Bowen 


EW (A) 

-^peakZ-^con 

FWHM (km s" 1 ) 


6.26±0.03 
1.3±Q.l 
919±11 


4.06±0.07 
1.3±Q.l 
99^16* 



For a single narrow Bowen line only. 



discrepancy, we used a weighted average of both ephemerides 
(assuming all quoted errors are correct) to derive the or- 
bital phase for our spectra. We show the results in TableQ] 
In order to estimate the phase uncertainty, we use a com- 
bination of the weighted formal error of two ephemerides 
(0.05 orbital phase) plus a systematic error of 0.09 orbital 
phase derived from the difference between the formal and 
observed error between the 2 ephemerides. This gives a to- 
tal minimum phase uncertainty on our spectra of 0.1 orbital 
phase. 

In order to determine the radial velocity for each spec- 
trum we fitted Gaussians to the two narrow NIII lines 
plus narrow component of He II A4686 simultaneously, where 
the rest-wavelength of each Gaussian and their separation 
were fixed. We then determined the off-set from their rest- 
wavelength using a minimum \ 2 technique. We also deter- 
mined the off-set using each line individually, but the dif- 




4620 4640 4660 4680 4700 4720 
Wavelength (A) 

Figure 3. Close-up of the Hell A4686 and Bowen region for a 
single observation obtained on 27 May 2004 (bottom) and an 
average of all observations taken on 18 June 2004 (top). 



ferences are negligible. However, the presence of the narrow 
stationary lines could have an influence on our estimated 
off-sets, especially close to orbital phases 0.5 and 1.0 where 
they blend with the moving lines. 

Fig.|6]shows our radial velocity measurements as a func- 
tion of orbital phase. Note that due to observing conditions 
we have averaged the 4 spectra around orbital phase 0.9 into 
two bins to increase the signal to noise. For illustration pur- 
poses we also included the 18 June spectra (indicated with 
filled circles and labeled as '18 Jun'), where we used the nar- 
row, stationary, lines to determine the off-set. Fitting a sine 
curve to the measurements, with the orbital period taken as 
fixed, it is clear that the 18 June lines do not originate from 
the same source in the binary frame, and we therefore do not 
include them in our radial velocity fits. From a fit through 
the resulting points we derived a velocity semi-amplitude of 
if em =247±8 km s~\ a mean velocity of 7=30±10 km s" 1 
and phase zero of </>o=0.88±0.05. We determined the errors 
on the fit by scaling the errors on the individual radial ve- 
locities by a factor of ~4 such that a xt=l is achieved. Note 
that phase zero was a free parameter, and does not include 
any systematic uncertainty in the ephemeris used which is 
expected to be of order ~0.1. We do notice a small jump in 
the radial velocities around orbital phase 0.33, due to the 
narrow stationary lines that are starting to blend with the 
wing of the moving components, but excluding these points 
did not change the results of our fit. 

A close inspection of the broad component of Hell 
A4686 shows that it moves in anti-phase with the narrow 
component. However, depending on the orbital phase either 
the wings or the core is dominated by the narrow component. 
We have tried to subtract this narrow component, but were 
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Figure 4. Close-up of H/3 during the June 2004 observations. 

not able to do this reliably and could therefore not measure 
the radial velocities for the broad component from this line. 
Instead, we investigated the other prominent emission lines 
and noted that the only line where a narrow component is 
not the dominating feature was H/3. However, also for this 
line we have to be cautious since there is some contamination 
from Hell A4859 as shown in Fig.(8] We applied the double 
Gaussian technique to the wings of this line (Schneider & 
Young 1980). As we moved away from the line core, we noted 
that the velocity semi-amplitude was stable at i\i=68±5 km 
s _1 , the zero-phase was at </>o=0.35±0.05 and the mean ve- 
locity around 7=130±20 km s (while keeping the orbital 
period fixed again). These values appear to be consistent 
with those from the broad component in Hell A4686, and 
indeed show a 0.5 phase difference between the broad and 
narrow components. However, we do note that the differ- 
ence in mean velocity between the narrow and broad lines is 
consistent with the difference between H/3 and Hell A4859, 
suggesting that the latter is contributing significantly to the 
blend. Furthermore, we also note that formally the fit is 
not acceptable (x 2 =458 for 36 d.o.f.), showing that we must 
treat this result with caution. 



4 DONOR STAR DETECTION 

We have detected narrow N III fluorescence lines in the opti- 
cal spectrum of Aql X-l during outburst that are very similar 
to the donor star fluorescence features observed in other sys- 
tems. Having derived the radial velocity curve for these fea- 
tures, we will now discuss the likelihood of these tracing the 
irradiated donor in Aql X-l. These lines are a dominant fea- 
ture of the Bowen blend, as was also the case in GX339— 4, 
ScoX-1, X 1822-371, V801Ara and V926 Sco (Hynes et al. 
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Figure 5. Trail of the Bowcn and He II A4686 lines clearly show- 
ing the movement of the lines due to orbital motion. Note that 
we have not phase-binned the spectra due to the sparse sampling. 



2003a, Steeghs & Casares 2002, Casares et al. 2003+2006). 
In these X-ray binaries the narrow lines are thought to orig- 
inate on the surface of the companion star. This is nicely 
illustrated in the eclipsing LMXB X 1822-371. Here, the ori- 
gin of the narrow Bowen lines can unambiguously be claimed 
due to the excellent agreement between the photometric 
and spectroscopic ephemerides (Casares et al. 2003). For 
the non-eclipsing LMXBs one has to rely on ephemerides 
based on determining the minima from lightcurves that are 
quite variable and far from sinusoidal. Again, when the num- 
ber of cycles is large enough one also finds an excellent 
agreement with the spectroscopic ephemerides, as was the 
case for V801 Ara, and also explained the apparent differ- 
ence between the photometric and spectroscopic phasing in 
V926Sco (Casares et al. 2006). We therefore think that the 
most likely place of origin of the narrow lines observed in 
AqlX-1 is also on the surface of the companion star. 

Unfortunately, due to the rather large uncertainty in 
the phasing (0.1 orbital phase), we cannot rule out the pos- 
sibility that the narrow lines in AqlX-1 do not originate 
on the companion star. However, for the other two poten- 
tial regions where these lines could originate, the accretion 
stream or impact point, it is not clear if it could produce 
such narrow emission lines as observed in AqlX-1. Compar- 
ing the Doppler maps of several LMXBs for which a sig- 
nature of these regions is observed, namely EXO 0748—678, 
X 1822-371 and XTE J1118+480, shows that the stream is 
observed either in absorption or the region is too extended 
to produce such narrow lines (Pearson et al. 2006, Casares 
et al. 2003, Torres et al. 2002). Furthermore, also in cata- 
clysmic variables the accretion stream or impact point are 
commonly observed, but again they are very broad and do 
not tend to have sinusoidal radial velocity curves. We there- 
fore think that the surface of the companion star is the most 
likely place of origin for these NIII narrow lines, although 
an up-dated ephemeris and a Doppler map of the Bowen 
region are needed to confirm this. 



4.1 Curious Features 

Narrow Bowen components have now been observed in grow- 
ing list of X-ray binaries (e.g. Steeghs & Casares 2002, 
Casares et al. 2003, Casares et al. 2006), and we can now 
also add AqlX-1 to that list. Interestingly, the narrow lines 
in AqlX-1 show a strong similarity with those observed in 
GX339— 4, also an X-ray transient but with a black hole 
primary that was confirmed by the observation of the nar- 
row Bowen components (Hynes et al. 2003a). In GX339— 4 
the narrow components were of a comparable strength to 
those in AqlX-1 and also dominate the Bowen region, while 
there was also an indication that a narrow component was 
contaminating the He II emission (Hynes et al. 2003a) . 

From the Bowen/He II radial velocity curve we have es- 
timated a systemic velocity of 30±10 km s _1 , while H/3 gives 
a value of 130T20 km s . This is most likely due to the 
presence of Hell A4859, and would suggest that it is domi- 
nating the blend. Such contamination has been observed in 
other X-ray binaries such as XTEJ2123— 058 (Hynes et al. 
2001). However, this does mean that we will have to be very 
careful with associating this semi-amplitude with the orbital 
motion of the neutron star, since this could also be distorted 
due to the blending of different components. Such caution 
is further enhanced by the fact that observations of, for ex- 
ample, X 1822-371 by Cowley et al. (2003) have shown that 
measuring the semi-amplitude of the compact object from 
emission lines can only give a first order approximation. A 
full orbital phase coverage of a non-blended line that will not 
have a strong narrow component during a future outburst 
of AqlX-1 will be needed to estimate the K- velocity of the 
neutron star. 

The most curious feature of AqlX-1 is the change of 
the emission lines during the June 2004 observations. Such 
a dramatic change in the morphology of the emission lines 
has been observed before in the X-ray nova GRO J0422+32 
(Shrader et al. 1996). When the source hardened in X- 
rays the Hell A4686 and the Bowen lines were suddenly 
greatly reduced in intensity while Ha and H/3 became much 
stronger. This change was interpreted as the "turning on" 
of the hard X-ray component. Although, no hard X-ray/soft 
gamma ray observations were performed for AqlX-1 around 
the time of our observations in order to test this (Bird, pri- 
vate comm.), the soft X-ray lightcurve (Fig.QJ does indicate 
a large change in the accretion activity between the May 
and June observations. 

Even if such "turning on" of the hard X-ray component 
occurred in AqlX-1 it is still strange that we do not observe 
the donor star anymore in the June 2004 observations. Since 
these observations were performed during the final stages of 
the outburst, i.e. a few days before it went back into quies- 
cence (Maitra & Bailyn 2005), a large fraction of the inner 
accretion disk might already have disappeared. This could 
have stopped the production of UV photons that are needed 
in the Bowen fluorescence process (McClintock et al. 1975), 
and thereby terminating the narrow lines coming from the 
companion. Unfortunately, this scenario cannot be verified 
with the limited spectral information that is contained in 
the ASM data. Another suggestion could be that the open- 
ing angle of the accretion disk has increased between the 
different observations, and is screening the donor star at the 
end of the outburst (see also Sect. 4.2). 
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Figure 6. Best fit radial velocity curve of the narrow Bowen blend 
components and He II A4686 as derived from a multi-Gaussian fit. 
Phase zero corresponds to inferior conjunction of the donor star 
according to the ephemeris derived in the text (and illustrated as 
the dashed curve) .Note that the error on the ephemeris is at least 
0.1 orbital phase. Also indicated as filled circles and labeled '18 
Jun' are the radial velocities derived from the June 2004 spectra 
that have been excluded from the fit. We have shown the curve 
twice for clarity. 
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Figure 7. Best fit radial velocity curve derived from H/3 based 
on the double Gaussian method. We have shown the curve twice 
for clarity. 



The other curious feature of AqlX-1 is the narrow, sta- 
tionary lines that are present in both the May and June 
observations. These narrow lines do not show any change 
in radial velocity between the observed spectra, and have 
a mean velocity of more or less 30 km s~ , which is inter- 
estingly comparable to the systemic velocity derived from 
the lines that arise on the donor star. AqlX-1 is not the 
first source where such low velocity emission lines are ob- 
served. For example, in the dwarf novae IP Peg, SS Cyg 
and U Gem among others (see Unda-Sanzana et al. 2006 
for an overview) also show such low velocity emission lines 
(Steeghs et al. 1996, Unda-Sunzana et al. 2006). Also in the 
LMXB SerX-1 such unusually narrow and stationary lines 
were present (Hynes et al. 2003b). This makes AqlX-1 an- 
other LMXB to become a member of the group that shows 
these low velocity emission lines, and clearly indicates that 
these lines are a long-lived and common feature that are 
independent of the spectral state of the source. Currently 
there is no satisfactory explanation for these low velocity 
emission lines, but it is becoming apparent that these lines 
are present in many compact binaries. Therefore, a convinc- 
ing explanation should be sought. However, it is clear from 
AqlX-1 that they are not connected to the companion, con- 



trary to the narrow and moving lines observed in the May 
2004 spectra. 

4.2 System parameters 

If we assume that we have detected a signature of the com- 
panion star, what can we say, in combination with previous 
observations, about the components of AqlX-1? We know 
that the compact object is a neutron star due to the pres- 
ence of type I X-ray bursts (Koyama et al. 1981), leading to 
a mass of the compact object <3.1M0. Based on photom- 
etry and spectroscopy of the optical counterpart of AqlX-1 
(correcting for the interloper star at a distance of 0.46 arc- 
sec) Chevalier et al. (1999) determined the spectral type of 
the secondary as K6-M0. Although several constraints on 
the inclination exist, and even contradict each other (Shah- 
baz et al. 1997, 1998, Garcia et al. 1999, Welsh et al. 2000), 
most were determined before it was known that the opti- 
cal counterpart consisted of a blend of two stars (Chevalier 
et al. 1999). Shahbaz et al. (1997) determined an inclination 
of ~50° using the width of the Ha line, while Garcia et al. 
(1999) argued that the inclination must be ~70° based on 
the amplitude of the variations in the outburst lightcurve. 
Finally, Welsh et al. (2000) used the de-blended optical data 
to determine that the inclination is >36°, while the absence 
of dips or eclipses gives an upper-limit on the inclination of 
< 70-80° (Paczynski 1971). 

From Fig.[U] we have determined a K em velocity of 
247±8 km s~ . However, since the narrow lines are expected 
to arise on the surface of the donor star, this value will 
be a lower-limit on the centre of mass velocity of the sec- 
ondary. Without using any other constraint on AqlX-1 we 
can determine that the mass function f(Mi)=M\ sin 3 + 
g) 2 >1.23±O.12M0 , only slightly smaller than the canonical 
value of a neutron star (IAMq). We note that an inclination 
^73°, which would be consistent with all other observations, 
pushes the mass of the neutron star in AqlX-1 above this 
canonical value. 

Following Casares et al. (2006) we decided to use the 
width of the emission lines to estimate a lower limit to the 
rotational broadening, Kotsini, of the donor star. If we as- 
sume a synchronously rotating companion, this gives an esti- 
mate of q and K 2 via Kotsini=0.4627 i 2g 1/3 (1 + q) 2/3 (Wade 
& Home 1988). In order to estimate Kotsini we created an 
average spectrum in the rest frame of the narrow compo- 
nents, and the Hell A4686 and H/3 regions are shown in 
Fig. [8] From this spectrum we estimated the FWHM of the 
NIII A4634.12 line at 147.8±36.7 km s" 1 (we did not use 
NIII A4640.64 since it is blended with NIII A4641.85). We 
also estimated the FWHM of NIII A4634.12 in all individ- 
ual spectra and find an average of 135.7±14.8 km s . We 
do note that there was no detectable systematic change in 
FWHM as a function of orbital phase. However, these values 
do include the effect of the intrinsic broadening due to the 
instrumental resolution of 70 km s . To take this into ac- 
count we broadened a strong line in our arc spectrum using a 
Gray rotational spectrum (Gray 1992) until we reached the 
observed FWHM. We assumed no limb-darking since the flu- 
orescence lines occur in optically thin conditions. We found 
that a rotational broadening of 96±15 km s _1 is required. 

We have to be very careful with our estimate for the 
rotational broadening. First of all, we assumed that all the 
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Figure 8. Top: Average spectrum of AqlX-1, around the Hell 
A4686 and Bowen region, in the rest-frame of the donor star. 
Bottom: Average spectrum of AqlX-1, around H/3, in the rest- 
frame of the donor star. 



emission of the narrow component is coming from the com- 
panion. However, we do note a variation in line strength 
(see Fig.[5j| and there might also be other broadening mecha- 
nisms, as was most likely the case in GX 339—4 (Hynes et al. 
2003a). In this source the widths of emission lines were too 
high (270-350 km s" 1 ) to be due to rotational broadening 
only, and we cannot exclude that some other broadening 
mechanisms are also present in AqlX-1. Furthermore, we 



also should take into account the effect of smearing due to 
the orbital motion on the emission lines. We estimated that 
it should be at maximum ~20 km s _1 for AqlX-1, and there- 
fore smeared our arc-line by this amount. We noticed that 
this changed our estimates on the rotational broadening by 
SCO. 6 km s _1 , small enough to be discarded. Finally, since 
the rotational broadening is close to the spectral resolution 
(70 km s -1 ), the optimal V ro tsim found could be very sen- 
sitive to the assumed template. We therefore also used the 
[OI] A557.4 sky-line of all the individual spectra to deter- 
mine the rotational broadening in this way. We notice that 
the variation in seeing has a small effect on the estimate 
for the broadening, but it was always ^2 km s _1 . Although 
the difference between the sky-line and arc-line is negligible, 
high resolution spectra of a ~full orbit will be needed to un- 
ambiguously determine the rotational broadening. Bearing 
in mind these caveats, let us investigate the implications of 
our Kotsini estimates on the neutron star mass. 

Combining Kotsini of 96±15 km s _1 with the so-called 
Jf-correction (that allows us to estimate the centre of mass 
velocity of the secondary using K eln ) we can determine both 
q and K% (Muhoz-D arias et al. 2005). However, we do note 
that this i^-correction is (weakly) dependent on the inclina- 
tion of the system, but strongly dependent on the opening 
angle of the accretion disk, both of which are unknown. We 
used the fourth order polynomial fits given by Mufioz-Darias 
(2005) to determine the /^-correction (and also the corrected 
mass ratio q) that gives the lowest estimate for the neutron 
star mass. We did notice that this is a monotonically de- 
creasing function with opening angle, i.e. the largest open- 
ing angle gives the smallest Tf-correction and therefore the 
smallest neutron star mass. On the other hand, the largest 
disk opening angle does give the largest mass ratio q. Assum- 
ing that 7^2=247±8 km s -1 , we estimate that <j<0.58 at 95% 
confidence (a disk opening angle of ~18°) corresponding to 
Afisin 3 i^2.2±0.3Mo , larger than the canonical neutron star 
mass. However, we do stress again that high-resolution spec- 
tra are needed to verify this result, but also note that thus 
far we have assumed no i\-correction (and thus underesti- 
mated K2) or taken into account the inclination, both of 
which will only increase the mass of the neutron star. 

Let us in addition speculate that the velocity derived 
from Fig.[7](68±5 km s~ x ) gives a first order estimate for the 
orbital motion of the neutron star. In this case (^=0.28) is 
close to that determined above (if K eul is close to the centre 
of mass velocity K2), and leads to a (maximum) disk open- 
ing angle of a~16°. Since AqlX-1 does not show eclipses, we 
can follow Paczynski (1974, 1983) to estimate an upper-limit 
on the inclination using cosi= 0.46(l+q)~ 1/3 . This leads to 
i<65° and together with q=0.28 gives an estimate on the 
neutron star mass of Mi>2.7Mq, suggesting a rather mas- 
sive neutron star. Unfortunately, since we cannot trust the 
derived radial velocity of the primary we cannot unambigu- 
ously conclude this, but if the K\ velocity could be confirmed 
and the inclination of the system constrained, it would then 
place strong limits on the equation of state of neutron star 
material (e.g. Lattimer & Prakash 2001). 

While estimating the ^-correction for the different disk 
opening angles, we did notice that an angle of ~14-15° gives 
a minimum mass for the neutron star >3TM©. This sug- 
gests that already for opening angles slightly smaller than 
the maximum opening angle derived above, the compact ob- 
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ject in AqlX-1 will exceed the maximum mass for a neutron 
star. Since we know that the compact object is a neutron 
star, this could mean that the opening angle during the May 
observations is close to the allowed maximum. If the open- 
ing angle of the accretion disk was slightly larger during 
the time of our June observations (either due to an asym- 
metry in the accretion disk or a real increase), this would 
completely screen the donor star and no narrow components 
would be observed anymore, explaining their absence during 
the observations at the end of the outburst. 



5 CONCLUSION 

We conclude that we have detected narrow emission lines 
in the Bowen region and Hell A4686 during an outburst of 
AqlX-1, that we interpret as coming from irradiated side 
of the donor star. If correct, this allows us to constrain 
its orbital velocity for the first time. This leads to a hard 
mass function for the neutron star of /(Mi) = 1.23±O.12M0 , 
and for an inclination ^73° exceeds the mass of a canonical 
1.4Mq neutron star. If we can trust the rotational broad- 
ening that we determined, the width of the emission lines 
places a lower limit on the mass of the primary of ^1.6Mq 
(at 95% confidence). This suggests that AqlX-1 harbours a 
more massive neutron star, warranting future observations 
exploiting the Bowen fluorescence lines in order to improve 
the constraints on its binary parameters. 

Better phase-coverage of the entire outburst with high- 
resolution spectra together with an updated ephemeris 
would not only constrain the orbital velocity of the compact 
object, but also help in understanding how the structure of 
the accretion disk evolves, and further our understanding of 
the nature of the low-velocity lines. Given the difficulty of 
observing AqlX-1 in quiescence, this paper shows that high 
resolution optical spectroscopy could be a promising way 
forward. 
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